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ABSTRACT 


Transportation of fluids in the industry especially in the oil and gas sector is a common phenomenon. Thus, analyzing the flow 
characteristics in elbow ducts is essential for a number of engineering applications, but the limited availability of theories and 
models to accurately calculate the pressure drop at predetermined velocities is somewhat intricate. Star-CCM+ v5.04 Computational 
Fluid Dynamics package was used to analyze the pressure characteristics at 5m/s, 10m/s, 20m/s and 40m/s velocity in a non-linear 
duct (elbow). The results obtained were compared to the calculated ones from the theoretical analysis for the same elbow ducts. It 
was found that the differences in pressure obtained were mainly as a result of performing the theoretical calculations using one 
dimensional simplified energy equation while the Computational Fluid Dynamics simulation was performed in two dimensions which 
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correlate more with what happens in a real flow. The influence of turbulence and boundary layer conditions which was taken into 


account more accurately by the software also contributed to the errors. 
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1. INTRODUCTION 


In the process of transporting fluid from one point to another through ducts or pipes which may either be linear or non-linear 
depending on the piping system, pressure drop is oftentimes encountered due to resistance to flow. Pressure gain/loss can also be 
encountered in ducts as a result of change in elevation between the starting and ending point, and may result in minor or major 
pressure difference (across the duct) which relates to a number of factors such as; pressure gain as a result of fluid head added 
through pumps, effects of friction acting between the fluid and the wall of the pipe, loss in pressure as a result of the change in 
elevation of the fluid, friction loss as the fluid flows through bends, pipe fittings etc. [1]. 

Considering the aforementioned factors, industries such as oil and gas encounter problems in attempt to analyze pressure 
difference across the piping network. Ellis and Joubert [2] carried out an investigation on turbulent shear flow in curved pipes by 
comparing the results obtained from two curved rectangular flow to the flow in straight pipe. Crawford et al. [3] successfully used 
frictional effects to predict a single phase pressure drop for turbulent flow at 90° bend, with no available models to separate the 
flow. Furthermore, Speddling et al. [4] in their analysis developed empirical correlations and also determine the pressure drop in 
elbow bends for single phase laminar and turbulent flow. 

Computational Fluid Dynamics (CFDs) is a branch of Fluid Mechanics that applies the principles of numerical analysis and data 
structures in analyzing pressure distributions in fluid flow [5]. Briley et al. [6] in their analyses used Navier-Stokes equation to 
numerical compute a turbulent flow in 90° square duct channels and pipe bends. The said differential pressure can also be 
determine by theoretical calculations, the use of gauges etc. but a more specialised tool such as CFDs may also be required when 
determining the fluid flow at non-linear 90° elbow bends [1, 4, 6,7] which is the primary focus of this study. 


2. RESEARCH METHODOLOGY 


Mesh convergence studies was carried out for different velocities but at 20m/s velocity, the solution was found to converge at a 
mesh size of 15mm as shown in Table 1. This mesh size was therefore used to run the analysis for other velocities. Surface mesh was 
used to mesh the surface of the model. Surface mesh re-triangulates an existing surface in order to improve the overall quality of 
the surface and optimize it for the volume mesh models. The Trimmer option was selected for the volume meshing. The mesh 
convergence study is presented in Table 1. The option provides a robust and efficient method of producing a high quality grid for 
both simple and complex mesh generation problems and produces predominantly hexahedral mesh with minimal cell skewness. 
Turbulence intensity is a measure of the intensity of turbulence of fluid flow. In flows within a closed conduit, the turbulence 
intensity varies from 1% to 10% [8]. 


Table 1 Mesh Convergence study 


Mesh Sizes Velocity 
(mm) (m/s) 
60 60 
50 50 
40 40 
30 30 
15 20 
10 10 


For the simulation, turbulence intensity was set at 10%. Turbulent length scale is a physical quantity describing the size of the 
large energy-containing eddies in a turbulent flow. In pipe flows, turbulent length scale can be estimated from the hydraulic 
diameter and has been shown to vary between 3.8% and 7%. For this simulation, the turbulent length scale was taken as 7% of 
hydraulic diameter. In the same manner, 5% turbulent velocity scale of the free steam velocity was used in the simulation. In this 
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study, theoretical calculation and star CCM simulation was used to determine the static and total pressure in non-linear flow system 


for 5m/s, 10m/s, 20m/s and 40m/s velocities and both the theoretical calculations and simulated results were compared for possible 
errors and difference. Mesh convergence study for the elbow pipe is presented in Table 1. The channel length of the pipe was 2.6m 
with pipe diameter of 16.8mm. Mesh control was applied to the model in order to establish some salient factors such as the element 
shape, mid-side node placement and element size. Curvature based mesh was also applied in order to refine all regions of higher 
curvature, as these regions are prone to stress formations due to the non-linear shape of the duct. These information are basically 
for the model development process, and can affect the accuracy of the model and subsequent analysis. 


2.1. Boundary Conditions 
Table 2 shows the boundary conditions obtained after the analysis for the respective velocities. 


Table 2 Boundary conditions 


Inlet Boundary Mass Flow 
; Mach Turbulence 
Velocity Layer rate(Inlet) : : 
number ; intensity 
(m/s) Thickness (kg/s) 
5 0.0128 0.06670 2.9 0.29881 
10 0.0186 0.06041 2.9 0.29881 
20 0.0312 0.05471 2.9 0.29881 
40 0.121 0.04955 2.9 0.29881 


This represents the total pressure (contour), Velocity (vector) and Static Pressure (Contours) for 5m/s, 10m/s, 20m/s and 40m/s 
velocity at corresponding mass flow rate of 2.9kg/s, 3.6kg/s, 5.4kg/s and 7.8 kg/s respectively. Flow output across the internal 
gradient of the 90° elbow duct for each velocity is presented in Figure 1 to 12. The turbulent length scale was 6% of the hydraulic 
diameter while turbulent velocity scale was 8% of the free steam velocity. 


2.2. Total Pressure (Contours) 


Total Pressure (Pad 
ooe2 


a OE43)4 3 73IIS 12.28) 16 SSS 20.839 
Figure 1 Plot of total pressure at 5m/s 

g Total Pressure (Fa) 

ss SSISE 20.594 35689 S170 66.353 sié13 


Figure 2 Plot of total pressure at 10m/s 
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f" Toto! Pressure (Pa) 
160.60 213.83 


ar 54.154 107.38 267.05 320.27 


Figure 3 Plot of total pressure at 20m/s 


bx 317.7) 5G2 &7 
Figure 4 Plot of total pressure at 40m/s 


2.3. Velocity (Vector) 


Velocity (m/s) 
24 0.67870 1.8070 2.9354 40637 §.1920 63206 


Figure 5 Plot of Velocity at 5m/s 
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a4 3.2777 5.1370 


Figure 6 Plot of Velocity at 10m/s 


JS 


f Velocity (m/s) 
Zé 0.068973 5.0847 10.100 15.116 20.132 25.147 
Figure 7 Plot of Velocity at 20m/s 
f Velocity (m/s) 
2s 16.265 23.056 29.845 35.637 43,427 50.218 


Figure 8 Plot of Velocity at 40m/s 
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2.4. Static Pressure (Contours) 


S 


Static Pressure (Fa) 
as 9.0835 5.3227 1.5657 21993 5.9502 97212 


JS 


| Static Pressure (Po} 
4 -59.164 -23.71T -8.2561 71947 22.647 36. 5G6 


Figure 9 Plot of static pressure at 5m/s 


Figure 10 Plot of static pressure at 10m/s 


Static Pressure (Po) 


2 133 84,58) -23 568 32.853 91370 159.29 


Figure 11 Plot of static pressure at 20m/s 
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as 


Static Pressure Pa} 
a # 531.32 -350.77 -120.12 Mast 347.42 571.73 


—v, 


Figure 12 Plot of static pressure at 40m/s 


2.5. Theoretical Analysis 

Analyzing the flow characteristics in air ducts is essential for a number of engineering applications, but the limited availability of 
theories and models to accurately determine the pressure condition at predetermined velocities is somewhat intricate. This section 
involves the application of existing equations in calculating pressure conditions in the elbow duct at different scenarios. 


2.6. Calculation of Static Pressure 
To find the inlet static pressure by theoretical calculation, Equation (1) was used; 


Po = P (inlet static) + 0.5(pV?) i 


where, 

Py = Total inlet pressure (Pa), 

Piinlet static. = Inlet static pressure (Pa), 

p = Density of air (kg/m?), V = Inlet velocity (m/s). 


Total Inlet pressure = Inlet static pressure + Inlet dynamic pressure (2) 
Hence; 
Pantet static) — Po = 05 (pV?) (3) 


Substituting the values obtained by theoretical calculations for total inlet pressure into Equation (1), the following values were 
obtained as shown in Table 3; 


Table 3 Values Obtained from Theoretical Calculations 


Inlet Velocity (V) (m/s) 


Density (p) (kg/m?) 


The total inlet 
pressure (P.) (Pa) 


Inlet static pressure 
Piinlet static) (Pa) 


5 1.2 56.900 41.9 
10 12 118.350 58.35 
20 1.2 356.724 116.728 
40 1.2 1285.956 325.956 


The Energy equation is given by Equation (4) 


2.7. Pressure loss obtained from 1D Formula 
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Py 
pg 


P2 


LE tps 4 a (4) 
29 al loss ~ ng 29 2 


where, 

P; = The total pressure at the inlet, 

P2 = The total pressure at the outlet which is equal the atmospheric pressure, 

V; = the inlet velocity, V2 = the outlet velocity, 

Z, = the distance between inlet and the datum point, 

Zz = the distance between outlet and the datum point, 

p = the density of air at 20°C 

For continuity, 

Qi = Q2 (5) 


where, 


Q; = mass flow rate at the inlet, 
Q> = mass flow rate at the outlet; 


Hence, the mass flow rate at the inlet is given by; 
PAV, = pA, (6) 


When the inlet and outlet area are the same, and the density are the same, then the inlet velocity is equal to the outlet velocity as 
shown in Equation (7); 


V, = V2 (7) 


Moreover, the outlet pressure is equal to the atmospheric pressure, which is equal zero as shown in Equation (8); 


we 
il 
Oo 


(8) 


By substitute the equation 7 and 8 into equation 4, Equation (9) was obtained as follow; 


P, P2 | (V2-V;)? 

og 7 (22 — 41) + Rigas t ag (9) 
Py 

bg = (Z2 —Z4) + Nioss(totat) (10) 
Py =pg X|[Z2— 21) + hiossctotan| (11) 


2.8. Pressure Losses at the Duct 
Three losses happen in the duct as the flow passed from inlet to outlet; these losses are major loss in horizontal length and vertical 
length, and minor loss in two bends. The following steps were adopted in the pressure loss theoretical calculations for 5m/s velocity. 


hictotat) = Aimajor loss x-direction) + Aymajor loss x-direction) + Ay(minor loss at bend) (1 2) 


pVD _ 12x5x04 


Re= = 
m 1.8 x10 


(13) 
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Re = 133,333.33 = 1.33 x 105 


The flow is turbulent since Re is greater than 4,000. The following represents the theoretical background of the pressure analysis. 
Major losses due to the horizontal length of the pipe. The head loss due to friction can be obtained; 


L v2 
Head loss (major loss x-direction) =fcre) x = x 2g (14) 


Major losses due to the vertical length of the pipe. The head loss due to friction is shown in Equation (15). 


Lo. V2 
Head loss (major loss Y-direction) = fre) x a x 2g (15) 


Minor losses due to the bends of the pipe. The head loss due to friction was obtained as follow: 


v2 
Head loss (minor loss at bend) =K, X 2g (16) 
Since there is two bends at 90° of the system 


y2 
Head loss (minor loss at bend) = 2 X K, X 2g (17) 


The mass flow rate was given by Equation 18; 
Q = pAV (18) 


where, 


Q = The mass flow rate (m?/s), 

A = The area for 1m width and 0.4 thick (m2), 
V = The velocity (m/s), 

p= The density of air (kg/m?). 


0.4m was used as the area of the duct since the simulation was done in 2D. The following steps were used in the pressure loss; 
theoretical calculations for 10m/s, 20m/s and 40m/s velocity. Turbulence kinetic energy for Standard k-epsilon turbulence model is 
given by Equation 19, while turbulent dissipation for Standard k-epsilon turbulence model is given by Equation 20; 


@ fou — (ug He) *] =p (5, 2 — 8.220) pe 2 (4, 44 ok) OM 
at (pk) + Ox; Jew (u 7 *) =| H, (si dx;  oh,tp 2.) Pe 3 (u, Ox; a pk) Ox; (19) 
0 (nee 0 ( i He) O€ 
a ag AP aoe, 
du; 2 du; du; 2 ; 10 du; 
Cor | Sy got — 2 (H, ga! + ok) Sy — Coop = — Cosp Ep, S22 + Coupe oe (20) 


p, k, €, u, 6, Si represents density, turbulent kinetic energy, turbulent dissipation rate, absolute velocity, turbulence Prandtl number 
and average tension rate, while Ce1, Ce2, Cez, Cegstandard k-epsilon Turbulence Coefficients. 


3. RESULTS AND DISCUSSION 


The following results shown in Table 4-9 were obtained from both the theoretical calculations and simulation of flow conditions in 
the 90° elbow duct, whereas, Figure 13-16 shows the graphical representation for various calculations. 
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Table 4 Results Obtained from Theoretical Calculations 


S/N Velocity (m/s) Reynolds number fire) u maior ee hi maior oa 
direction) (m) direction) (m) 
1 5 133,333.33 0.01688 0.4840 0.0753 
2 10 266,666.67 0.01477 1.6933 0.2634 
3 20 533,333.33 0.01305 5.9845 0.9309 
4 40 1,066,666.67 0.01164 21.3656 3.3235 
Table 5 Total Inlet Pressure and Mass Flow Rate from Theoretical Calculations 
S/N Velocity hi (minor loss at hy (total) (Mm) Z2-Z\ Total Inlet Mass flow rate 
(m/s) bend) (m) (m) pressure (Pa) (kg/s) 
1 5 1.2742 1.8335 3 56.900 24 
2 10 5.0968 7.0535 3 118.350 48 
3 20 20.3874 27.3028 3 356.724 9.6 
4 40 81.5494 106.2385 3 1285.956 19.2 
Table 6 Data from the CFDs Simulation 
Inlet Velocity cules Inlet Pressure Outlet Inlet mass Flow | Outlet mass Flow rate 
(m/s) pane (pa) Pressure (pa) rate (kg/s) (kg/s) 
5 5.328105 5.94 0 24 2.4 
10 10.54154 21.91 0 48 48 
20 20.97109 81.58 0 9.6 9.6 
40 41.65581 285.66 0 1.92 19.92 
Table 7 Values Obtained from the CFDs Simulation 
Inlet Velocity | Inlet total pressure Outlet total Inlet static Outlet static 
(m/s) (pa) pressure (pa) pressure (pa) pressure (pa) 
5 20.94 17.03 5.94 0 
10 81.92 66.67 21.91 0 
20 321.58 263.87 81.58 0 
40 1245.66 1041.12 285.66 0 


Table 8 Total Inlet Pressure (Theoretical Calculation) and Inlet Pressure (Simulation) 


Velocity (m/s) Total Inlet pressure (Pa) Total Inlet pressure (Pa) % Error Between theoretical 
theoretical calculation Simulation calculation and Simulation 

5 56.90 20.94 63.20 

10 118.35 81.91 30.79 

20 356.72 321.58 9.85 

40 1285.96 1245.66 3.13 


Table 9 Comparing both Total Inlet Pressure for Theoretical Calculation and Simulation 


; Static Inlet pressure from Static Inlet pressure from 
Inlet Velocity (V) (m/s) j : ; j Percentage of Error 
theoretical calculation simulation 
5 41.9 5.94 85.82 
10 58.35 21.91 62.45 
20 116.728 81.58 30.11 
40 325.956 285.66 12.36 
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Total Inlet Pressures vs Velocity 


Inlet pressure (pa) 


5 10 20 40 
Hand calculation 56.9 118.35 356.724 1285.956 
Simulation 15mm Mesh size 20.94 81.91 321.58 1245.66 


Velocity m/s 


Hand calculation —— Simulation 15mm Mesh size 


Figure 13 Total Inlet pressure (Theoretical Calculations and Simulation) versus Velocity 


Error (difference between the total Inlet pressures of simulation and theoretical 
calculations) vs Velocity 
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Inlet pressure for Hand calculation vs Inlet pressure for Simulation at 15mm Mesh size 


Figure 14 Difference in Total Inlet pressure (Theoretical Calculation and Simulation) versus Velocity 
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Static Inlet Pressure vs Velocity 


NY WH W 
— 
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150 


Total Inlet pressure (pa) 


Velocity m/s 


Hand calculation —— Simulation 15mm Mesh size 


Figure 15 Static Inlet Pressure (Theoretical Calculation and Simulation) versus Velocity 


Percentage Error (difference between the Static inlet pressures of simulation anad theoretical 
calculations) vs Velocity 


100.00 
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85.82 


Percentage of error 


Velocity 


Inlet pressure for Hand calculation vs Inlet pressure for Simulation at 15mm Mesh size 


Figure 16 Difference between the Static Inlet Pressures versus Velocity 


Figure 13 shows total pressure against inlet velocity for the results of theoretical calculations and simulations, while Figure 14 
shows the plot of error between the total inlet pressures and the inlet velocities. Figure 15 depicts the results of plot of static inlet 
pressures for simulations and theoretical calculations against inlet velocity while Figure 16 shows the plot of error between the static 
pressures and inlet velocity. As shown in Figure 13 and 15, it can be observed that the pressure increased with increase in velocity. 
The increase is exponential because the velocity term in the pressure equation is squared. As shown in Figure 14 and 16, it can be 
observed that the error decreases as velocity is increased from 5m/s to 40m/s. The difference between the results obtained from 
theoretical calculations and that obtained from CFDs analysis can be explained to be as a result of performing the theoretical 
calculations in one dimension while the simulation was done in two dimensions. 

The one dimension analysis simplifies flows assuming velocity and fluid energy to be the same throughout the flow. The two 
dimensions is closer and more accurate representation of what is happening in reality where the fluid velocity varies from point to 
point. Also the effect of other flow parameters such as the boundary layers where not taken into account in one dimension 
calculation, while its effect was taken into account in the computational fluid dynamics analysis. In general, for simple geometries 
such as the one used in this analysis, the one dimension pressure loss formula can be used to get a general idea about flow 
conditions while more detailed analysis can be done in two dimensions or three dimensions for more accurate and detailed results. 
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Analysis in computational fluid dynamics is affected by the number of grid points (cells) generated to solve the computation. The 


number of cell generated is a function of the mesh size. 

Generally, as the number of cell increases, the results obtained become more accurate while the computational time increases 
also. However, as the mesh size is made finer and the number of cells increased, a point is reached when the results obtained is not 
or is marginally affected by the mesh size. At this point, the mesh is said to have converged. The results obtained at this point are 
usually taken as the solution of the computation. Due to the length of time required to obtain solutions using fine mesh sizes, initial 
analyses were done using coarse mesh. The mesh size was gradually refined until convergence was achieved. Another important 
parameter which affected the results obtained from the simulation is the number of iteration to convergence. The iteration steps 
were increased until the results obtained were stabilized (i.e., the results were no longer changing with time). In general, flow in a 
two dimensional plane is considered as a special case of a three dimensional if the geometry is symmetrical in one coordinate [9]. 
However, experiments have shown that two dimensional models give a very close approximation to three dimensional models for 
symmetrical model [10, 11]. 


4. CONCLUSION 


Three dimensional flow dynamics in non-linear duct (90° elbow) at predetermined velocities were analyzed using Star CCM+, and 
the flow parameters and characteristics were successfully determined. The results obtained from the simulations were compared 
with the ones obtained from theoretical calculations. The error between the two results were explained to be as a result of 
performing the theoretical calculations using one dimensional energy equation, whereas, the simulations were carried out in two 
dimensions which was a better representation of what happens in real life situation. The effects of boundary layer conditions which 
was taken into account in the simulations was also a contributing factor to the errors obtained. 
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Nomenclature 


CFDs Computational Fluid Dynamics 
Po Total inlet pressure 
Pressure 
p Density 
V Velocity 
Zi Distance between inlet and the datum point 
Z2 Distance between outlet and the datum point 
Q Mass flow rate 
A Area 
Kg/s Kilogram per seconds 
m/s Meters per seconds 
Re Reynold number 
Ky Frictional Factor 
K Turbulent kinetic energy 
€ Turbulent dissipation rate 
Absolute velocity 
fo} Turbulence Prandtl number 
Si Average tension rate 
hy Head loss 
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